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The possibility for our visible world to be a 3-brane embedded in a multidimensional bulk is at 
the heart of many theoretical edifices in high-energy physics. Probing the braneworld hypothesis 
is thus a major experimental challenge. Following recent theoretical works showing that matter 
swapping between braneworlds can occur, we propose a neutron-shining-through-a-wall experiment. 

We first show that an intense neutron source such as a nuclear reactor core can induce a hidden 
neutron flux in an adjacent hidden braneworld. We then describe how a low-background detector 
can detect neutrons arising from the hidden world and quantify the expected sensitivity to the 
swapping probability. As a proof of concept, a constraint is derived from previous experiments. 

PACS numbers: 11.25.Wx, 12.60.-i, 28.20.-v 


I. INTRODUCTION 

As suggested by several extensions of the standard 
models of particle physics and cosmology, our world could 
form a three-dimensional space sheet (a brane) embedded 
in a larger bulk with more dimensions As for most 

theories beyond the Standard Model of particle physics, 
new effects predicted by the braneworld hypothesis at 
both the high-energy and the precision frontiers can be 
investigated. Particle colliders, such as the Large Hadron 
Collider, attempt to excite new degrees of freedom at 
high energies |4j, |5J - like the Kaluza-Klein excitations 
of particles in the extra dimensions - whereas precision 
experiments at low energies attempt to find tiny signals 
induced by the new physics [M3- For instance, in the 
context of braneworld scenarios, the compactified extra 
dimensions would modify the inverse-square law of grav¬ 
ity at short distance. A great variety of experimental 
techniques have been developed to search for a possible 
modification of gravity from subatomic to macroscopic 
distance scales (see Ref. jb] for a recent review). 

In the present work, we are interested in some pecu¬ 
liar low-energy effects induced by the existence of other 
branes in the bulk. In recent theoretical works [T3 - [T7j | . 
it has been argued that usual matter could leap from 
our braneworld to a hidden one, and vice versa (see 
Fig. 1). This matter swapping between two neighbor¬ 
ing braneworlds could be triggered by magnetic vector 
potentials, either of astrophysical origin or generated ex¬ 
perimentally. In particular, we have considered neutrons 
oscillating between a state where they sit in our brane 
and a state corresponding to a neutron located in the 
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other brane. The matter-swapping probability would os¬ 
cillate at high frequency r] and small amplitude p. More¬ 
over, neutrons are well-known versatile tools which are 
already used to test other concepts such as axion [l8j or 
mirror particles [l9l - l23j . 

In a previous work, we have shown that such oscil¬ 
lations could have affected experiments measuring the 
neutron lifetime. The existing experiments set a con¬ 
straint on anomalous neutron disappearance at the level 
of p < 7 x 10 -6 [13 • It is amusing to note that a 
small tension has very recently appeared between the 
neutron lifetimes measured with two different methods 
[24l |. which could be interpreted as anomalous neutron 
disappearance. In addition, such a constraint on neutron 
disappearance could also be useful to constrain some ap¬ 
proaches related to the big bang nucleosynthesis [13 ■ 

In the present paper, we propose an experiment to in¬ 
vestigate matter swapping between branes by looking at 
the appearance of neutrons from a neighboring brane. 
The concept is similar to light-shining-through-a-wall ex¬ 
periments (see Ref. Q for a recent review of this topic) 
where photons from an intense light source would con¬ 
vert into a sterile state (dark photons or axion-like par¬ 
ticles) which could pass through a wall, then would con¬ 
vert back to photons. In a neutron-shining-through-a- 
wall experiment, a very bright source of neutrons and a 
low-background neutron detector separated by a wall are 
necessary. Neutrons would swap into a sterile state - the 
state where they are located in another brane - which 
would be free to cross the wall. Then, the reappearance 
of neutrons into the detector situated behind the wall 
is checked. Although they are in principle very similar, 
there is an essential difference between light and neu¬ 
trons shining through walls. As far as light is concerned, 
the oscillation frequency is supposed to be rather slow 
and the conversion of photons into a sterile state builds 
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up coherently over long distance, while the oscillation is 
rather fast for neutrons and their conversion into sterile 
states results from the successive collisions of neutrons at 
nuclei. 

In section II, the theoretical aspects and the phe¬ 
nomenology of the neutron dynamics in a two-brane uni¬ 
verse at low energy are reviewed. The conditions lead¬ 
ing to matter swapping between branes are given. The 
strength of the ambient magnetic vector potential (which 
drives the matter exchange between branes) is discussed, 
as well as the environmental conditions that could pre¬ 
clude the swapping to occur. In section III, we describe 
how the neutron diffusion in the reactor moderator cre¬ 
ates a neutron flux in another brane and how this hidden 
neutron flux can be detected. Finally, the sensitivity of 
the suggested experiment is estimated, and as a proof of 
the concept, a constraint is derived from previous exper¬ 
iments. 


tion [T3-IT~H : 

! 4(t) =(H ° +H "" ) (t)' (1) 

where the indices ± allow one to discriminate the two 
branes. and ip- are usual Pauli spinors correspond¬ 
ing to the wave functions in the (+) and (—) branes, 
respectively. The unperturbed Hamiltonian Ho is block 
diagonal and describes the usual nonrelativistic evolution 
of two uncoupled spin-1/2 fields in the magnetic fields B± 
and the gravitational potentials V±, in each brane: 

H ° = ( H 0 + H_ ) ■ < 2 > 

U ± = ~A- Iin a-B ± + V ± . (3) 


II. PHENOMENOLOGY OF BRANE MATTER 
SWAPPING 


A. Dynamics of swapping 



S=l/g 


The nondiagonal part H cm of the Hamiltonian describes 
the coupling between the branes, 


H 


cm — dfj'ri 


0 

ia- {A + — A_} 


-ia- {A + — A_} 
0 


(4) 


where A+ and A_ correspond to the magnetic vector 
potentials in the branes (+) and (—), respectively. The 
same convention is applied to the magnetic fields B± 
and to the gravitational potentials V±. g, n is the mag¬ 
netic moment of the neutron. H cm implies that matter 
exchange between branes depends on the magnetic mo¬ 
ment and on the difference between the local (i.e., on 
a brane) values of the magnetic vector potentials, g is 
the coupling strength between the matter fields of each 
brane. 

In the following section, the case of an ambient mag¬ 
netic potential with an astrophysical origin is considered 
(section |HB]). Let A amb = A ambj+ - A amb _ be the dif¬ 
ference between the ambient magnetic potentials of each 
brane. Assuming that g, n B± -C V±. i.e., the magnetic 
fields in the branes can be neglected (in particular as¬ 
suming that V x A arnb = 0), then by solving the Pauli 
equation, the probability for a neutron initially localized 
in our brane to be found in the other brane is fl 2 j 


FIG. 1. Naive view of the neutron swapping between two 
branes. A neutron n in our brane can be transferred into 
a hidden brane (neutron n ) under the influence of a suit¬ 
able magnetic vector potential A. g represents the coupling 
strength between the branes related to the effective distance 
S between each brane. 

In previous works [M3 , it was shown that, in a uni¬ 
verse containing two parallel braneworlds, the quantum 
dynamics of a spin— 1/2 fermion in the presence of an 
electromagnetic field can be described by a two-brane 
Pauli equation at low energies. In electromagnetic and 
gravitational fields, the dynamics of a free neutron, i.e., 
outside of a nucleus, is described by the following equa- 


p = Sin2 ((V2) , (5) 

where g = \V + — V-\/K and H = ggi n A arnb /A P is the 
instantaneous probability for matter swapping between 
branes. Equation ([5]) shows that the neutron in the po¬ 
tential A amb undergoes Rabi-like oscillations between the 
branes. Note that the swapping probability is indepen¬ 
dent of the neutron spin direction [3] . As detailed in 
previous papers |14H16I|. the environmental interactions 
(related to V±) are usually strong enough and the oscil¬ 
lations are suppressed by the factor fl /77 << 1 . As a 
consequence, in a nucleus, a neutron is fully frozen in its 
brane due to the large nuclear potential. 
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B. Ambient magnetic vector potential 

The overall ambient astrophysical magnetic vector po¬ 
tential A amb was previously assessed in the literature 
[HI]. A amb is the sum of all of the magnetic vector po¬ 
tential contributions related to the magnetic fields of as¬ 
trophysical objects (planets, stars, galaxies, etc.) since 
B(r) = V x A(r). As a rule of thumb, A ss DB where 
D is the distance from the astrophysical source and B 
is the typical field induced by the object. At large dis¬ 
tances from sources (for instance, close to the Earth), 
A amb is almost uniform (i.e. V x A nm b ss 0) and can¬ 
not be canceled with magnetic shields [2^]. Now, Eqs. 
© and © show the dependence of the swapping effect 
against A amb = A amb , + - A amb _, i.e. the difference 
between the vector potentials of the two braneworlds. 
Since A„ m j _ depends on unknown sources in the hidden 
brane, we cannot assess its value. Then, A amb should be 
considered as an unknown parameter of the model. Nev¬ 
ertheless, the expected order of magnitude of A amb can 
be roughly constrained by A ambt+ in our visible world 
(since A amb results from a vectorial difference, it seems 
quite unlikely that A amb can fortuitously fall to zero). 
Galactic magnetic field variations on local scales towards 
the Milky Way core ( A amb ss 2 x 10 9 T m) are usually as¬ 
sumed [26]. By contrast, the Earth’s magnetic field leads 
to 200 T m while the Sun contributes 10 T m [2(|. By 
contrast, intergalactic contributions were expected to be 
about 10 12 T m [Hi! ( see also our previous papers for a 
more detailed discussion [Hl-Qjj]). Anyway, for now A amb 
is a parameter fairly bounded between 10 9 T m and 10 12 
T m. 


C. Environmental potential 

If one considers free neutrons shielded from magnetic 
fields, only gravitational contributions are relevant. Be¬ 
cause r] = \V+-V-\/h, it is difficult to assess the value of 
gH as it results from a scalar difference involving the un¬ 
known gravitational contribution V- of the hidden world. 
Therefore, g appears as an effective unknown parameter 
of the model and could reach weak values of a few eV up 
to large values around 1 keV. Indeed, estimations given 
in previous works [H-[l3] suggest that V + could be of the 
order of 500 eV due to the Milky Way core gravitational 
influence on neutrons. Note that the Sun, the Earth, 
and the Moon provide lower contributions of about 9 eV, 
0.65 eV, and 0.1 meV, respectively. At last, it must be 
emphasized that g is time dependent due the motion of 
the Earth around the Sun at the lab scale. Between the 
Earth’s aphelion and perihelion, the gravitational energy 
(due to the Sun) of a neutron varies from 9.12 to 9.43 eV. 
Of course, a time dependence could have many different 
origins. For instance, the particle motion relative to an 
unknown mass distribution in the hidden brane should 
be considered. However, it is unlikely that the Earth 
is ’’close” enough to a hidden mass distribution that is 


large enough to induce a significant energy time depen¬ 
dence on a time scale of about one year or one day. The 
time dependence is then mainly induced by the Earth’s 
motion around the Sun, with a relative variation in g 
of about 0.3%o in one year. Such a variation could be 
detected through an annual modulation of the swapping 
probability (see section ITlIDD . 

D. Neutrons as a sensitive probe 

With high-energy particle colliders, the braneworld hy¬ 
pothesis can be investigated if the brane thickness is in 
the range £ « 10 “ 19 m corresponding to the TeV scale 
( hc/t ; ~ 1 TeV). Colliders are blind to the Planck scale 
10 -35 m. By contrast, experiments at lower energies us¬ 
ing high-intensity neutron sources could reveal a multi- 
brane world through effects induced by the interbrane 
coupling g , which can be approximated by jhi ] 

g oc ( 1/0 exp {-kd/C) , ( 6 ) 

where d is the real distance between each brane in the 
bulk and k is a constant of the model [l2j] . Now, as 
an illustration, let us consider (for instance) a coupling 
constant g «10 ~ 3 m _1 . Such a value is consistent with 
present experimental bounds E3- For branes at the TeV 
scale, the above value of g is reached for d ~ 50£. Now, 
if one considers branes at the Planck scale, the coupling 
constant remains unchanged for d « 87£. As a conse¬ 
quence, while brane physics could be invisible for collid¬ 
ers, it could be observed in low-energy experiments using 
neutrons if a second brane exists close enough to ours. 

At last, we note that neutrons are more suitable than 
electrons, protons, or atoms for such a purpose phi ]. In¬ 
deed, for a charged particle, the Hamiltonian © also 
contains the usual terms, 

H p ± = — — A± • P + |A±| 2 , (7) 

m Zm 

with P = —ih\7. This implies that the term Hg in Eq. 
(151) is then supplemented by large terms proportional to 
qA amb which strongly freeze the oscillations. Consider¬ 
ing, for instance, A amb xs 200 T m (see section III Bl) 
and a proton with a kinetic energy of about 10 eV, the 
first term in the Hamiltonian © contributes 9 MeV to 
Hi 7 . For such values, the amplitude of the oscillations 
is suppressed by 8 orders of magnitude compared to the 
neutron case. Now, let us consider atoms. Though they 
are neutral, atoms are endowed with an instantaneous 
electric dipole moment (IEDM) d. Obviously, accord¬ 
ing to the time average, we must verify that (d) = 0 
and (d 2 ) 7 ^ 0 since IEDM results from quantum fluctu¬ 
ations of atomic orbitals. The London dispersion forces 
between atoms result from interactions between these in¬ 
stantaneous dipoles. Then, the Hamiltonian ([3]) must be 
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supplemented by terms which derive from Eq. 0: 

H<i,± - -A± ■ (8) 

where the fluctuation time is about At ~ K/Ei , where 
Ei is the ionization energy of the atom (Ei ~ 10 eV). 
Provided that At is larger than the period T of the Rabi 
oscillation of an atom between two branes, hr] in Eq. 
m is then supplemented by a huge term proportional to 
(l/At)A am (, • d. Then, for d about 10 ~ 30 C-m (0.3 D) 
(a typical value for atoms) and A am b ~ 200 T m only 
(for instance), we still get hr] ~ 20 MeV. With these 
values, the oscillations are strongly damped by 9 orders 
of magnitude compared to the neutron case. As a result, 
the neutron is a good candidate to test matter swapping 
between branes, since it is devoid of global charge or any 
electric dipole moment. 

E. Collision-induced neutron swapping probability 

As a consequence of the environmental potential i]h, 
the neutron oscillations present weak amplitude and high 
angular frequency of the order i]/2. Due to the fast oscil¬ 
lating behavior, one can approximate the swapping prob¬ 
ability by its time-averaged value p = (P) (see Eq. (0), 
such that 


When freely propagating, the neutron can be described 
as a superposition of two states: a neutron in our brane 
vs a neutron in the other brane. When colliding with 
a nucleus situated in our brane, the interaction acts as 
a measurement and the neutron collapses either in our 
brane with a probability 1 — p or in the other invisible 
brane with a probability p. In the following sections, the 
swapping probability p is considered as the relevant mea¬ 
surable parameter. By contrast, bounds on the coupling 
parameter g depend on the knowledge of galactic mag¬ 
netic potential fields and on the ambient gravitational 
fields. 


III. A 

NEUTRON-SHINING-THROUGH-A-WALL 
EXPERIMENT 

As explained above, in a two-brane universe, neutrons 
have a nonzero probability to escape from our brane into 
another brane at each collision. Therefore, a nuclear 
reactor, where the neutron density is very high, would be 
a very intense source of hidden neutrons. These neutrons 
could escape the reactor and be detected, after having 
swapped back to our brane, with a standard neutron 
detector located near the reactor (see Fig. 2). In the 
present section, the expected magnitude of the hidden 



Heavy water tank 


FIG. 2. Sketch of the experimental device. The source 
of possible hidden neutrons is a nuclear reactor core (here 
for instance, the Institut Laue-Langevin facility in Grenoble, 
France). The neutron detector using helium-3 gas is consid¬ 
ered to detect neutrons which could emerge from a hidden 
world. The detector is embedded in a shield to reach a very 
low background. 

neutron flux as a function of the swapping probability is 
discussed and the sensitivity which can be reached with 
such an experiment is estimated. For sake of clarity, the 
details of our calculations only appear in the Appendices. 


A. Induced neutron flux in the hidden world 

Let us assume that our Universe is made of two mu¬ 
tually invisible braneworlds and let us also consider a 
neutron flux <f> + inside a nuclear core. Our aim is to de¬ 
termine the intensity of the hidden neutron flux (in 
the vacuum of the hidden braneworld) in the core vicin¬ 
ity. For a given volume element of the reactor, since each 
neutron collision creates a hidden neutron with probabil¬ 
ity p, the hidden neutron source is proportional to the 
macroscopic elastic cross section T E of the reactor mod¬ 
erator and to the neutron flux More specifically, we 
get the source term corresponding to the number of gen¬ 
erated hidden neutrons per unit volume and unit time 
(see Appendix [All : 

S- = (10) 

where p is given by Eq. 0. Equation (1101) is derived 
by using the matrix density approach p7i. l28j| related to 
Eq. 0 as shown in Appendix [Al From this source term, 
we deduce the hidden neutron flux at the position r 
by considering the solid angle and integrating over the 
reactor volume V : 

$_(r) = -p / -- 2 T E (r')<f>+(r')d 3 r'. (11) 

ott Jv |r - r | 

The relation m shows that any signal should decay 
as 1/D 2 when the distance D between the reactor and 
the detector increases. This is an important issue to 
discriminate a hidden neutron signal from the common 
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neutron background in the reactor vicinity. Indeed, such 
a background mainly depends on local secondary sources. 


B. Principle of hidden neutron flux detection 

We now want to build a detector that is able to mea¬ 
sure the hidden neutron flux. Neutron detectors are 
based on the detection of charged particles emitted af¬ 
ter neutron absorption. In the present case, we suggest 
the use of a gaseous detector, such as the usual helium-3 
or boron trifluoride (BF 3 ) neutron detectors. The mecha¬ 
nism of hidden neutron capture can be described through 
the approach of Feinberg and Weinberg [27[, which was 
also used by Demidov, Gorbunov, and Tokareva [28j to 
describe the positronium oscillations in a volume full of 
gas in the mirror matter concept. Details of the calcu¬ 
lation are shown in Appendix [EJ We can then compute 
the event rate T detected in our brane against the swap¬ 
ping probability p (see Appendix iBl). For monochromatic 
neutrons, we get the intuitive but not obvious result 

r = ip£ A $_P, (12) 

where V is the volume of the detector and £U is the 
macroscopic absorption cross section related to BF 3 or 
He-3. For a continuous energy spectrum, the event rate 
is obtained by integrating Eq. ED over the spectrum. 
The main difficulty comes from the background which is 
usually high in such an environment. The detector must 
be shielded from neutrons. 


C. Proposal at the ILL reactor 

For our experiment, we propose to use the Institute 
Laue-Langevin (ILL) reactor (see Fig. 2) where the core 
(diameter of 40 cm and height of 80 cm) is surrounded 
by a heavy water tank (diameter of 2.5 m and height 
of 2 m) which will be our hidden neutron source. It 
can be shown that the major contribution to the hid¬ 
den neutron flux comes from the heavy water tank, while 
neutrons are strongly absorbed by light water and other 
surrounding materials. The thermal neutron flux inside 
the heavy water tank is modeled by using a point-like 
source <F + ( r ) = 5'exp(— r/L)/(4irDr) } where L = 116 
cm and D = 0.57 cm are the length and coefficient of dif¬ 
fusion for heavy water, respectively. S is fitted to match 
with the known neutron flux in the heavy water tank 
[Ill] with S = 6.0 x 10 17 neutrons/s. For heavy water, 
the macroscopic elastic cross section is He = 0.4 cm -1 . 
The present rough model is reliable enough to discuss 
the experimental concept. Indeed, it is sufficient to as¬ 
sess the correct magnitude of the induced hidden neutron 
flux. Of course, any discussion of the later experimental 


results will require a detailed computation of the neu¬ 
tron flux in each location inside of the reactor, as well as 
a consideration of the various materials surrounding the 
core. 

We plan to use a cylindrical helium-3 detector (see Fig. 
2) with a volume of 36 cm 3 and a gas pressure of 4 atm. 
This detector could be located at 10 meters from the 
center of the core. The detector will be shielded by a 
47 r-box made of a (at least) 3-cm thick borated rubber 
(40% boron content) to capture thermal neutrons, sup¬ 
plemented by a polyethylene cover with a tunable thick¬ 
ness to moderate epithermal and fast neutrons. Indeed, 
hidden neutrons should provide a characteristic constant 
counting rate which does not depend on the shielding 
thickness. 

Using Eqs. ED and ED, the estimated rates are 
shown in Fig. 3 (black solid line). Considering the previ¬ 
ous constraint p < 7 x 10 - 6 M (vertical blue solid line), 
we see that the expected event rate could be as high as 
36 kHz, which is easily detectable even without shielding. 
Actually, we think we can reach an upper rate of about 
10 mHz or even 1 mHz f30j - [32| with suitable shields and 
detector (horizontal red solid line). We can then expect 
to reach a new upper constraint for the swapping proba¬ 
bility p of about 10 -9 . Such a constraint would be better 
than the constraint p < 7 x 10 “ 6 0 by at least 3 orders 
of magnitude. 



Swapping Probability p 


FIG. 3. (Color online). Black solid line: Expected event rate 
T against the swapping probability p. Horizontal red solid 
line: Expected threshold of the background noise (1 mHz). 
Vertical blue solid line: Known upper limit of the constraint 
on the swapping probability p < 7 x 10 " 6 0. Black dashed 
line: Expected event rate E against the swapping probability 
p by considering the experiment in Ref. j30(. Horizontal red 
dashed line: Background noise (0.3 mHz) in the experiment 
in Ref. [30). Vertical blue dashed line: Threshold of the 
rough constraint on the swapping probability (p < 5 x 10 -9 ) 
considering the experiment in Ref. [3Cj]. 

In Fig. 4, as an example, we show the resulting ex¬ 
pected bound on the coupling constant g between our 
brane and an invisible one. The values of g are given 
against the gravitational constraint 77 from Eq. © • The 
existing constraint is given (shaded domain above the 
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short-dashed blue line) and is compared to the expected 
results of the neutron-shining-through-a-wall experiment 
(black solid line) by assuming A am b « 2 x 10 9 T m [TtI . 
For the best rate constraint, we improve the constraint on 
g by 2 orders of magnitude. By contrast, a constraint on 
the detected signal lower than 10 events per second (see 
Fig. 3) still allows p < 10~ 7 , i.e., we improve the con¬ 
straint on g by 1 order of magnitude. As a consequence, 
the present experiment can easily improve the constraint 
on the coupling between our world and invisible ones. 



FIG. 4. (Color online). Expected experimental limits for 
the coupling constant g against the confinement energy for 
some constraints on the swapping probability p for A arn i, = 
2 x 10 9 T m. Red-grey domains are excluded, p = 7 x 10 -6 
corresponds to our previous best experimental constraint, p = 
5 x ICC 9 corresponds to the expected constraint deduced from 
the experiment in Ref. [30l i|. 


D. Yearly time-dependent drift of the swapping 
probability 

Let us now briefly underline the possible consequence 
of the potential gH time evolution due to the Earth’s 
revolution around the Sun, as described in section III Cl 
From Eq. © p varies in time as A p/p = 2Ag/g. Now, 
from Eqs. m and m r varies as Ar/r = 2A p/p. 
Then, we deduce that the event rate T varies in time as 
Ar/r = AAg/g, i.e., AT/T « 10~ 3 over six months if 
we consider the values from section III Cl If it is expected 
to detect a time-dependent drift of the swapping prob¬ 
ability, for a given duration At of the experiment, the 
number of detected neutrons N = TAf and its instru¬ 
mental uncertainty AN ~ y/N imply that enough neu¬ 
trons must be detected to allow A N/N to be lower than 
AAg/g. In the present case, for month-by-month mea¬ 
surements, we should have N > 10 6 detected neutrons 
per month at least. Assuming At ~ 20 days per month, 
this leads to T « 1 Hz. As a consequence, considering the 
above-mentioned detector, a time-dependent drift of the 
swapping probability p can be measured provided that 
p> 3.7x10-®. 


E. A constraint as a proof of concept 

A constraint can be suggested from previous experi¬ 
ments j30j-[32j. Indeed, in different contexts, these exper¬ 
iments required low-noise neutron detectors to detect ul¬ 
tracold neutrons generated by conventional sources. But 
since these detectors were in the nuclear reactor neighbor¬ 
hood, they should have been sensitive to hidden neutrons 
as well. As a consequence, considering the recorded back¬ 
ground of these experiments, we can already derive a first 
rough constraint on the swapping probability, and thus 
on the coupling constant between two adjacent branes. 
We consider the favorable conditions introduced in Ref. 
pS0| . The detector volume was 500 cm 3 with a partial 
helium-3 pressure of 15 mbar. The distance between the 
detector and the nuclear core was 16.5 m. The back¬ 
ground was 0.3 mHz [3C|. The resulting constraint is 
shown in Fig. 3. Using the previously introduced reactor 
model, the expected hidden neutron flux F against the 
swapping probability p is shown by the black dashed line. 
The neutron background is shown by the horizontal red 
dashed line. As a result, this leads to a new constraint 
for the swapping probability p such that (see vertical blue 
dashed line in Fig. 3) 

p<5x 10” 9 . (13) 

Such a value can be related to a constraint on the cou¬ 
pling constant between our visible braneworld and a hy¬ 
pothetical invisible one. This constraint is given by the 
red dashed line in Fig. 4. Obviously, these constraints 
need to be confirmed or improved by a dedicated experi¬ 
ment such as that described in the previous sections. As 
a consequence, the present constraint on p does not pre¬ 
clude the measurement of a yearly time-dependent drift 
of the swapping probability. Nevertheless, it shows the 
feasibility of the experiment introduced in the present 
work. 


IV. CONCLUSION 

We have proposed an experiment in which neutrons 
can be used to test the braneworld hypothesis. Using 
an intense flux of neutrons behind a wall and utilizing 
a properly shielded neutron detector, reappearing neu¬ 
trons could be detected as a proof of the existence of 
hidden braneworlds. Three typical signatures can be 
examined to verify the reliability of a measured signal 
as a real exotic phenomenon. First, for thick enough 
shielding, the number of detected neutrons must not de¬ 
pend on the shield thickness. Second, when the distance 
to the core increases, the number of detected neutrons 
must decay as an inverse-square law against the distance 
between the nuclear core and the detector. At last, a 
yearly modulated time-dependent drift of the number of 
detected neutrons could be observed for a not too low 
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swapping probability. By contrast, without direct de¬ 
tection, such an experiment at least allows one to con¬ 
strain the existence of braneworlds. The experimental 
constraint on the matter-swapping probability between 
branes could be improved by 3 orders of magnitude with 
the proposed experiment. Then, a rough constraint on 
the swapping probability was proposed to show the fea¬ 
sibility of a neutron-shining-through-a-wall experiment. 


ACKNOWLEDGMENTS 


where V g ^ and V, h - are the gravitational potentials 
felt by the neutron in each brane. We assume that 
<7 ■ A arn b = AambO'z and we restrain ourselves to a spin- 
up state without loss of generality (the matter-swapping 
probability does not depend on the direction of the mag¬ 
netic vector potential or the spin state flill - ITfil ]). Then we 
can write 


H = 


~£i A + V 9,+ ~ inh ^ 

-£. A + Vg _)' 


(A4) 


The authors thank Karin Derochette for critical read¬ 
ing of the manuscript. 


Appendix A: Shielded nuclear core as a hidden 
neutron source 

In a previous work El , it was shown that ultracold 
neutrons colliding with a solid wall can escape to another 
hidden brane. In this situation, since ultracold neutrons 
present a wavelength roughly larger than the typical dis¬ 
tance between atoms of the wall, the wall appears as a 
continuous medium for neutrons [333]. Such a medium 
can be described through the Fermi potential [33] • Here, 
we want to estimate the production rate of hidden neu¬ 
trons due to collisions between fast (or thermal) neutrons 
and materials surrounding the nuclear core, i.e. mainly 
the heavy and light water tanks. In this context, since 
neutrons cannot be considered in the ’’optical” domain, 
the media cannot be considered as continuous, and we 
must consider collisions between neutrons and individ¬ 
ual nuclei. Then, our aim is to find a relation between 
the visible neutron flux and the hidden neutron flux 
<!>_ induced by the collisions between visible neutrons 
and nuclei in our brane. More specifically, we look for a 
source term, 


with Qh = gp n A amb . The behavior of the neutron flux 
is described through the use of the matrix density p [27[ , 

-p = -ih~ 1 ('Hp-p'H' l )+I c , (A5) 

where I c is the collisional integral which describes the 
collisions between neutrons and nuclei in the medium, 
with 


I c = nv J F( 0 )pFl( 0 )dn. (A 6 ) 

v is the mean relative velocity between neutrons and 
molecules and n is the density of molecules (nuclei) in 
matter. We have 


m = 


(m o 
loo 


(A7) 


which describes the scattering of neutrons by the 
molecules and where 0 is a scattering angle. The sec¬ 
ond diagonal term is equal to zero since we assume there 
is no matter in the hidden brane. 

We also define Ti. =H+C, with 


C = 


—2nnvtif(0)/k 0 \ 

0 0 ) ’ 


(AS) 


5_=A'$ + , (Al) 

where the constant K in the source term will be assessed 
in the following. It will be shown that K = (l/2)pE j g, 
where T,e is the macroscopic elastic cross section of the 
materials of the reactor and >I> + is the neutron flux in 
the whole reactor. To do this, considering the matrix 
density approach, we follow the Feinberg-Weinberg ap¬ 
proach [27|, which was also used by Demidov, Gorbunov, 
and Tokareva [28[ to describe other kinds of oscillations 
in a volume occupied by matter. 

Neutrons are described through the two-brane Pauli 
equation and we neglect the neutron decay process: 

ihd t ^ = Hf, (A2) 

such that 

_ f 2 m A 4“ Vg,+ idpi-iCb ' A. arnb 

y * A. amb T Vg — 


(where k is the neutron wave vector), which accounts for 
the presence of matter in the system. Let us use the 
definition 

p=( P+ . r + lS ). (A9) 

y r — is p- J v ' 

We also define 

TL'l) TIV 

wr = 4n— Re(/(0)), and wi = 47 t—I m(/(0)) (A10) 
k k 

with the optical theorem related to the cross section cr tot : 

<Jtot = cr E + cr i = -y-Im(/(0)) = —, (All) 

k nv 

where cje and o/ are the elastic and inelastic (absorption) 
cross sections, and 

<TE = J \mf dQ. (A12) 



We then obtain the following system from Eq. (IA5I) : 

—p+ = —2f2r — nvajp + , (A13) 

^P- = 2f2r, (A14) 

-^r = (77 - (1/2 )w R ) s—(l/2)rw I +£l (p+ - p_), (A15) 

= - (g - {1/2)w r ) r - (l/2)sw I: (A16) 

with hg = Vg t + — V g -. In the following, we consider 
a statistical set of neutrons initially localized in our 
visible brane and emerging from the nuclear core, i.e., 
p + (t = 0) = 1 and p_(t = 0) = 0. If we consider Eq. 
(IA13D . we obviously deduce that the neutron population 
in our brane must mainly decrease due to absorption, i.e., 
neutron leakage into a hidden brane must be very weak 
to be consistent with known physics. So, we must verify 
that 2f Ir -C nvt 7 ip+. Now, if (for instance) we consider 
thermal neutrons in heavy water, and if we consider 12 
with A am b = 2 x 10 9 T m and g = 10 -3 m _1 , we get 
nvai ~ 7.3 s _1 and 12 ss 2 x 10 14 s _1 , i.e., 2f2 nvai. 
As a consequence, the previous condition must imply that 
r <C p+- Since neutron leakage must be weak, we also 
deduce that p- <C p+- At last, assuming that g wr 
and using Eq. (IA11I) . Eqs. (IA13D - (IA16I) can be recast in 
the following convenient form: 

^ P+ = —nvcr I p+, (A17) 

-37 P- = 2f2r, (A18) 

at 

—r = r]s-(l/2)nv<j t otr + flp+, (A19) 

-778 = —ryr - (l/2)?w<r tot s. (A20) 

at 

At this stage of the problem, it is relevant to underline 
that from Eqs. (IA19I) and (IA20I) we can easily deduce 

= — (l/2)nva tot m + f2cosdp + , (A21) 

where we have set r + is = me lS . Without coupling, i.e., 
if 12 = 0, we deduce from Eq. (IA21I) that the coherence 
terms of the density matrix (i.e., its off-diagonal terms r 
and s) must exponentially decay as exp(—(1/2 )nv<Jtott)- 
As a consequence, the collisional dynamics should sup¬ 
press the quantum coherence precluding the neutron 
swapping between branes. By contrast, Eq. (IA21I) as well 
as Eqs. (IA19I) and (IA20I) also show us that the coupling 


12 can prevent decoherence. This allows stationary coher¬ 
ences such that dr/dt = ds/dt = 0, with r and s different 
from zero. The steady state is quickly achieved due to 
the above-mentioned exponential behavior. Let us call X 
the distance covered by the neutron in the medium. The 
length X cannot be equated to the crossed thickness L in 
the medium. Since neutrons are localized in our visible 
brane, they can be diffused by molecules or the nuclei of 
atoms. X is the sum of the lengths of the various straight 
lines covered by the neutron, and thus L = X aj/3 <je- 
Now, for instance, since na to t ~ 0.4 cm -1 for heavy wa¬ 
ter, the stationary coherences occur at X = 4 cm, i.e., 
L = 0.2 mm from the neutron source, a distance which 
must be compared to the 105-cm thick heavy water slab 
in the ILL’s reactor. By contrast, nai = 4 x 10 -5 cm -1 
in D 2 O, i.e., slowly varies in the heavy water region 
(see Eq. (IA17II L Then, in the stationary coherences hy¬ 
pothesis, Eqs. (IA19I) and (IA20D become a linear system 
of two equations with two unknowns (r and s) which is 
trivially solved. We get 

_1_ ( ( l/2)nva t otXLp+ \ 

(1/4) (nvatot) 2 + g 2 V ~ r I^P+ J 

(A22) 

From Eqs. (IA18I) and (IA22D we deduce 

^ P- ~ (1/2 )pnva to tP+, (A23) 

where we have set p = 2Q 2 /g 2 to be consistent with Eq. 
©• Let us now consider an initial local neutron flux 
4>o = uov where uq is the initial neutron density. The 
neutron fluxes in each brane are given by 4>-|- = 4>op±- 
Then, from Eq. (IA23I) we deduce 

<9 t <F_ = (1/2 )pnv<7 t ot®+- (A24) 

Since the neutron flux $_ = u-V, where U- is the local 
neutron density in the hidden brane, then 

d t u- = (1/2 )pna tot $+. (A25) 

Since U- is now local, Eq. (IA25I) must be supplemented 
by a divergence term V • j_ to account for the local be¬ 
havior of the neutron current j_ = zt_v. We then deduce 
the continuity equation for neutrons in the second brane, 



V • j_ + d t u- = (l/2)p?rcr tot $+, (A26) 

which is the continuity equation endowed with a source 
term S -, 


S- = (1/2 )pna tot A>+. (A27) 

Obviously, for heavy water, atot. ~ &E and we retrieve 
Eq. (EH) since tkje = E e■ Now, as a striking result, 
we note that neutrons propagating in the hidden brane 
are not absorbed since we expect vacuum rules in the 
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other world. As a consequence, while visible neutrons 
have been stopped by the water surrounding the nuclear 
core, hidden neutrons can propagate freely far away from 
the reactor. Then we can completely suppress the visible 
neutron flux in our brane while keeping a constant flux of 
hidden neutrons in the invisible brane. A hidden neutron 
source from a nuclear core surrounded by a relevant shield 
can then be built. 


Appendix B: Hidden neutron detector 


We consider a gas stored in a vessel as a neutron de¬ 
tector. More specifically, we consider a neutron flux in 
the second brane and we expect to detect neutrons which 
arise from the second brane in our braneworld. The dis¬ 
tance between the atoms in the gas is much larger than 
the neutrons’ wavelength. As a consequence, neutrons 
see a set of scattering objects instead of a continuous 
medium. We follow the same approach as before, and 
start with Eqs. (IA13I) to (IA16I) . but now we consider a 
statistical set of neutrons initially localized in the invis¬ 
ible brane, such as it can be written p+(t = 0 ) = 0 and 
P-{t = 0) = 1. 

Looking at Eq. (IA14I) , we see that the neutron popula¬ 
tion in the hidden brane mainly decreases due to neutron 
leakage into a hidden brane. Following a similar hypoth¬ 
esis as in the previous section, we can consider that 2 f Ir 
must be very weak and we assume that dtp- ~ 0 , i.e., 
is almost constant at the detector scale. We can also 
assume that <C p~ and ij 3> wr. Then (using Eq. 
(IA1 ID j . Eqs. (IA13I) - (IA16I) can be recast in the following 
convenient form: 


—p+ = -214r - nvcr[p +1 


(Bl) 


P- 


1 , 


(B2) 


—r = 77 s - (l/ 2 )rm;( 7 tot - 14p_, (B3) 


4 s = -pr - (l/2)snvatot- (B4) 

at 

Obviously, atot, cri, and cte are now the cross sections 
related to the gas and n is the density of gas molecules. 


Now, from Eqs. (1B3I) and (1B4I) we can easily deduce 
that 


—m = — (l/2)nvcrtotm — flcosdp-, (B5) 

where r + is = me 18 is set. Without coupling, i.e., if 
14 = 0, we deduce from Eq. (IB5I) that the coherence terms 
of the density matrix (i.e., its off-diagonal terms r and 
s) must exponentially decay as exp(— (l/2)nvatott)- As 
a consequence, the collisional dynamics should suppress 
the quantum coherence precluding the neutron swapping 
between branes. By contrast, Eq. (1B5I) as well as Eqs. 
(1B3I) and (1B4I) also show that the coupling 14 preserves 
from decoherence. Nevertheless, in contrast to the case 
described in the previous section, stationary coherences 
such that dr/dt = ds/dt = 0 (with r and s different from 
zero) cannot be achieved here due to the short length of 
the detector. Indeed, assuming for instance a He-3 detec¬ 
tor with a pressure of about 4 atm, we get ncr tot « 0.53 
cnr 1 for thermal neutrons. Now, the distance X = vt 
crossed by a neutron can be equated to the length L of 
the detector since we can assume that the hidden neu¬ 
trons are not significantly diffused by the nuclei of gas 
molecules or atoms in our brane. Let us consider (for in¬ 
stance) a detector length of about 10 cm. At half length, 
we get exp(— (l/2)nva to tt) ~ 27 %, i.e., the coherences 
are not negligible in at least half of the detector. As a 
consequence, we cannot apply the stationary coherences 
hypothesis here. 

Nevertheless, Eas. (IB3D and (IB4I) constitute a simple 
nonhomogeneous first-order differential equation that is 
easy to solve since p_ « 1, i.e., since is constant. Let 
us set 


( r \ p -(i/2)nv* tot t ( cos( 77 1 ) sin (rjt) \ 

\ s J y — sin(r)t) cos(r]t) J 

If we insert Eq. (IBGI) into Ecis. (IB3l) and (1B4D . 


(' 

(B 6 ) 
we obtain 


Ji/ 2 )nva tot t ( c.os(j]t) - sm(r)t) \ f -14 p_ \ 
dt\SJ \ sin (rjt) cos (??t) J ^ 0 J ’ 

(B7) 

from which we deduce that 



—14 


f* cos(7 ? ^ , )e ( 1 / 2 ) ™' fftott 'd4 , 
f* sin(r]t')e^ 2 '> nvatott 'dt’ 


(B 8 ) 


Using Eas. dBGI) and (IB 8 I) . Eq. (IBID can be rewritten 
as 


4 -P+ = 214V f e -(!/ 2 )nva tot (t-t') 

dt J o 

—nvaip+, (B9) 


which leads to 
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p + = 2 Q 2 p_ 


P- [ [ e - ( 1 / 2 ) ™ fft ° t(t '- t " ) e- n,KT/(t - t ' ) cos {p{t' - t”))dt''dt'. 

Jo Jo 


(BIO) 


From the properties of the density matrix p, we can de¬ 
duce the emerging neutron flux <E>„ which can be detected 
from the hidden current: = <F_p + . Considering He-3 

mixtures or BF 3 , we assume that the inelastic cross sec¬ 
tion ay stands for the neutron capture probability and 
the production of protons, which can be recorded. Then, 
the rate Y of events recorded per second by the detector 
in our brane is given by 


where Sd is the effective area of the detector and L is the 
length of the gas vessel, such that L/v is the time dur¬ 
ing which hidden neutrons travel in the detector volume. 
The length L can be equated to the length of the de¬ 
tector since we can assume that the hidden neutrons are 
not significantly diffused by the nuclei of gas molecules 
or atoms. In addition, absorption prevails on diffusion 
processes (ay <je in helium-3 or BF 3 ). 


rL/v r 

r= / nvai / <& v dSdt 

Jo Js d 

r L / v r 

= / nvaip + / &-dSdt, 

Jo Js d 


(Bll) 


Assuming that p ^ nvatot., and setting p = 2Yl 2 /p 2 and x 


vt, from Eqs. (IB 101) and (IB11I) we get 


T = ip s pJJ^JJJ- [ e~ naix [ [ e ( 1 / 2 ) niJtot ( x "~ X Je naix ' cos {{p/v) (x" — x')) dx"dx'dx (B12) 

v Jo Jo Jo 


with tp s = f s < F_(iS'. 


Due to the fast oscillating cosine term, and since p ^$> 
nvatot , Eq. (IB12I) reduces to 


detectors we simply get the intuitive but not so obvious 
result 


F - (1/2 )ip s p 



l){e~ n,JlL 


1 ) + na tot L 


(B13) 

For pure heliunr-3, relevant argon-helium-3 mixtures, or 
BF 3 , we verify that ay ay;, and then for long enough 


r ~ (l/2)j3noyV$_, (B14) 


where V = SdL is the volume of the detector and nay = 
E a is the macroscopic absorption cross section of the gas. 
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